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Electric current and spacial displacement due to trembling motion [Zitterbewegung (ZB)] of elec- 
trons in graphene in the presence of an external magnetic field are described. Contributions of both 
inequivalent K points in the Brillouin zone of graphene are considered. It is shown that, when the 
electrons are prepared in the form of wave packets, the presence of a quantizing magnetic field B 
has very important effects on ZB. (1) For B 7^ the ZB oscillations are permanent, for 5 = they 
are transient. (2) For B 7^ many ZB frequencies appear, for B = only one frequency is at work. 
(3) For B 7^ both interband and intraband (cyclotron) frequencies contribute to ZB, for B = Q 
there are no intraband frequencies. (4) Magnetic field intensity changes not only the ZB frequencies 
but the entire character of ZB spectrum. An emission of electromagnetic dipole radiation by the 
trembling electrons is proposed and described. It is argued that graphene in a magnetic field is a 
promising system for an experimental observation of Zitterbewegung. 

PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 
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I. INTRODUCTION 



The Zitterbewegung (ZB, trembling motion), first de- 
scribed by Schrodinger with the use of the Dirac equation 
for free relativistic electrons in a vacuum [l| , has been in 
recent years a subject of great theoretical interest since 
it has been shown that this phenomenon should exist in 
many systems in solids 0, i, i, S, i 0, S i, 0, [O, d, d • 
If one deals with two or more interacting energy branches, 
an interference of the upper and lower energy states gives 
rise to the ZB also in the absence of external fields. A 
formal similarity between two bands interacting via the 
k-p terms in a solid and the Dirac equation for relativistic 
electrons in a vacuum allows one to a pply similar theo- 
retical methods for both systems (see jl4l. Hsj). Most of 
the theoretical ZB treatments for semiconductors used as 
a starting point plane electron waves. However, Lock [l^ 
in his important paper observed that, since such a wave 
is not localized, it seems of a limited practicality to speak 
of rapid oscillations on the average position of a wave of 
infinite extent. Using the Dirac equation Lock demon- 
strated that, when an electron is represented by a wave 
packet, the ZB oscillations do not sustain their amplitude 
but become transient. The disappearance of oscillations 
at large times is guaranteed by the Riemann-Lebesgues 
theorem as long as the wave packet is a smoothly varying 
spatial function. The transient character of the trembling 
motion was demonstrated in our recent papers ^12,. ,17]. 
Since the ZB is by its nature not a stationary state but 
a dynamical phenomenon, it is natural to study it with 
the use of wave packets (see These have become 

of practical use with the emergence of the femto-second 
pulse technology jjj]. Recently, the transient trembling 
motion was proposed for ultra cold atoms 0, [2l|, for 



'Electronic address: [Tomasz. Rusin@centertel.pl| 



photons in two-dimensional photonic crystals |22l | and for 
Ramsey interferometry [2^ . Most recently, an actual ob- 
servation of an acoustic analogue of ZB in a macroscopic 
two-dimensional sonic crystal was reported [23 |. 

The trembling motion of charge carriers in solids has 
been described until present for no external potentials. 
On the other hand. Lock has remarked that, when an 
electron spectrum is discrete, the ZB oscillations do not 
have to disappear with time. In the present work we 
consider the trembling motion of electrons in solids in 
the presence of an external magnetic field, see [H, 
The magnetic field is known to cause no interband elec- 
tron transitions, so the essential features of ZB, which re- 
sults from an interference of positive and negative energy 
states of the system, are expected not to be destroyed. 
On the other hand, introduction of an external field pro- 
vides an important parameter affecting the ZB behavior. 
We show that, indeed, the ZB in a magnetic field is not 
damped in time. 

We consider the ZB phenomenon in monolayer 
graphene. This material has recently become an im- 
portant subject of investigation in the condensed matter 
physics and its applications [13, HI, |2§|. In particular, 
the charge carriers in graphene are considered to repre- 
sent massless Dirac fermions. In our approach we look 
for observable ZB phenomena. The situation we describe 
seems to be the most promising case for an experimental 
observation of the trembling motion in solids considered 
until present. 

The paper is organized as follows. First, we consider 
the Hamiltonian for electrons in graphene in a magnetic 
field, its eigenvalues and eigenfunctions. Second, we cal- 
culate carrier velocities and their averages taken over a 
Gaussian wave packet. Third, we give results for the 
ZB of electric current and electron position and empha- 
size the essential features introduced by the presence of a 
magnetic field. Finally, we propose and describe electro- 
magnetic radiation emitted by the trembling electrons. 
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The paper is concluded by a summary. In Appendices 
we discuss additional points related to the subject. 



II. PRELIMINARIES 

We consider a graphene monolayer in an external mag- 
netic field parallel to the z axis. As shown in both contin- 
uum and tight-binding models, there exist two inequiv- 
alent points K\ and K2 of the Brillouin zone (BZ) of 
graphene. The Hamiltonian for electrons and holes at 
the Kx point is (sO, El 











(1) 



where u « 1x10® cm/s is a characteristic velocity, tt = 
•p — qA is the generalized momentum, A is the vector 
potential and q is the electron charge. We remark that 
the above Hamiltonian is not symmetric in p^, and f)y. 
Using the Landau gauge we take A = {—By, 0,0), and 
for an electron q = — e with e > 0. Since in the above 
gauge H does not depend on x, we take the wave function 
in the form 



(2) 



Introducing the magnetic radius L — ^Jh/eB and the 
variable = y / L — k^L, we have 



H 



uh 







L \-i + d/d^ 



(3) 



Defining the standard raising and lowering operators for 
the harmonic oscillator a = -I- d/d£,)/^/2 and a+ — 
— d/d^)/^/2, the Hamiltonian becomes 



H = -hn 



a 
a+ 



(4) 



where the frequency is 

Next we determine the eigenstates and eigenenergies 
of the Hamiltonian H. Using a two-component function 
V' = (V'1,'02), we have 



E'ljj2. 



(5) 



Expressing ipi by -02 from the first equation we obtain 
from the second equation: h^Vi^a^ atp2 = The 
solution is ■!/'2 = |n), where \n) is the n-th state of the 
harmonic oscillator, and the energy is 



Ens — shrt^/n. 



(6) 



Here n — 0, 1, . . ., and s = ±1 for the conduction and 
valence bands, respectively. Formula ^ was experimen- 
tally confirmed for graphene in many works [28. . ,32, . .33.] . 
The functions (r|n) are given by 



(r|n) 



1 1 



(7) 



where C„ = y^2'"-n\y/TT , and H„ are the Hermite poly- 
nomials. From Eq. ([5]) we obtain 

Wi^ F — = r- — = -s|n-l). (8) 



E 



Each eigenstate |n) is labelled by three quantum num- 
bers: oscillator number n, energy branch s, and wave 
vector kx- The complete function is 



|n) = \nkxs) 



~s\n- 1) 



(9) 



For n = 0, the first component in Eq. ([9]) vanishes and 
the normalization coefficient is l/\/27r. 



III. VELOCITIES. ZITTERBEWEGUNG 

We want to calculate the velocity of charge carriers de- 
scribed by a wave packet. In order to do that we first cal- 
culate matrix elements (/|n) between an arbitrary two- 
component function / = (/",/') and eigenstates A 
straightforward manipulation gives 



{f\n) = -sF:_,+Fl 



(10) 



where 
in which 



g^{kx,y)e-i^"Rn{Ody, (11) 



P{x,y)e'^^-dx. (12) 



The superscript j — u,l stands for the upper and lower 
components of the function /. 

The Hamilton equations give the velocity components: 
Vi{0) = dH/diTi, with i = x,y. We want to calculate 
averages of the time-dependent velocity operators Vi(t) 
in the Heisenberg picture taken on the function /. The 
averages are 

m = ^e'^»'*/^/|n'>(n>,(0)|n)(n|f)e-'E°'/^ (13) 



where the energies and eigenstates are given in Eqs. ([6|) 
and The summation in Eq. (fT3|l goes over all the 
quantum numbers 



J j dkxdk'^Y^T.- 

nn' n fi ' a a' 



(14) 



We calculate a contribution to the velocity from the 
point Ki of the Brillouin zone. The matrix elements 
(n'|i;j,(0)|n) and (n'|t;a:(0)|n) can be shown to be 

{n'\Vy{0)\n) = —Sk^^k'{-sSn',n-l + s'Sn',n+l), (15) 
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(n>.(0)|n) 



:4,,fc' {sSn'^n-l + s'Sn'^n+i). (16) 



It is seen that the only non- vanishing matrix elements of 
the velocity components are those with the final states 



J 



n' = n ± 1. Putting the above matrix elements into Eq. 
P3p . we finally obtain after some manipulation for the 
Ki point of BZ 



Vyit) = w ^ y+ sin(w^t) + u ^ y„ sin(wf i) + ^ ^+ cos«i) + ^ A„ cos(wf i), 

n— n— n— 



oo 



v.{t) = uY,B+ cosiiu^J) + uY,B^ cos(wf t) + m ^ T+ sin(w^t) + ^ r„ sin(wf i), 

I 



n=0 



n=0 



(17) 
(18) 



where 



n=3 



^ n 

± _ 



T 

in which 



j-tU,U j-tU,U JJ^^^ TT^'^ 

_1_j-tU,U j-tU,U , JJ^^^ JJ^-^ 



Tjl,u _|_ Tr«,t 7-rl,'U 



Tju,l I 7-r(,u _|_ 7-r«,i _|_ jji.. 



jl.U 



(19) 



(20) 



The superscripts a, /3 refer to the upper and lower com- 
ponents, see Eqs. (fTU]) . The velocity averages must be 
real values. For example, if both g^{kx,y) {j — u,l) in 
Eq. are real then also both Fl{kx) given by Eq. pT|) 
are real, and there is = U^m- As a result the last 

two terms in Eqs. ([T7)) and ([TS]) vanish. 

The time dependent sine and cosine functions come 
from the exponential terms in Eq. (|13p . The frequen- 
cies in Eq s. ((Tf)) and (HH) are = ri(Vn + 1 - \Ai), 
w,f = 17(Vn -I- 1 -I- V^), where fl is given in Eq. (g]). The 
frequencies correspond to the intraband (cyclotron) 
energies while frequencies correspond to the inter- 
band energies, see Fig. 1. The intcrband frequencies are 
characteristic of the Zitterbewegung because the trem- 
bling motion is caused by an interference of states corre- 
sponding to the positive and negative energies [H, HH] . 
The intraband (cyclotron) energies are due to the band 
quantization by the magnetic field and they do not ap- 
pear in field- free situations (see (3. Isl. IgI. [T3| ) . 



IV. GAUSSIAN WAVE PACKET 

We take the function f{x, y) in the form of a Gaussian 
wave packet having an initial nonvanishing momentum 
POx = fikox 



f{x,y) = 



1 



^Tidxdy 



an 
ai 



(21) 




FIG. 1: The energy dispersion E{k) and the Landau levels 
for graphene in a magnetic field (schematically). Intraband 
(cyclotron) and interband energies for n' = n±l are indicated, 
see text. The basic energy is Ml = ^/2hu/ L. 



In the above model the upper and lower components of 
/ differ only by the coefficients a„ and a;, respectively. 



There is at 



9{kx,y) 



1. Then [see Eq. p2| ] 



-e 2- 



ai 



(22) 



This gives [see Eq. (fTTj) ] 



K{kx) 



'¥i(k^~ko.)\~hKD^ H„(fc,c) 



(23) 



4 



Here D = L"^ / J + dl,c = / ^ - d^, and 

n/2 



After some manipulations, we finally find (see [3a |) 



X (1 - (cQf) 




where Q = l/y^dl + D^, W = d.^DQ, and F = dlko^Q. 
Thus in case of a Gaussian wave packet we obtain the 
coefficients U^'^ in the form of analytical sums. For the 
special case dy — L, there is simply 



1/^,13 



where P = ^ -\- ^L"^ . In the above model the coeffi- 
cients J7"j'fj are real numbers, so that and in Eqs. 
(fT7| and (fT8| vanish. A sum rule for f/"'^ is given in 
Appendix A. 



V. RESULTS AND DISCUSSION 

In Fig. 2 we show calculated current components 
ii = —ev{t) as functions of time for different values of 
the initial wave vector k^x- The calculations were car- 
ried out for the Ki point of BZ taking = ai = l/v^ 
and a constant magnetic field of lOT. We assumed a cir- 
cular wave packet dx = dy = L — 83.13A, the coefficients 
^mn were calculated using Eq. ([2S]) . Figure 2a shows the 
results for ko^ = 0. It can been seen that, if there is no 
initial momentum, only jxit) is nonzero. The main fre- 
quency of oscillations is loq = f2, which can be interpreted 
either as luq — ri(-\/rn~T— y^) or ojq = ^l{^/n + 1 — ^/n) 
for n — 0. Frequency ujq belongs to the intraband (cy- 
clotron) set, while uj^ belongs to the interband set (see 
Fig. 1). Somewhat unequal amplitude of oscillations 
means that other frequencies also come into play, so that 
we already deal with the trembling motion. For fcox = Oi 
the asymmetry between Vx{t) and Vy{t) comes from the 
above mentioned asymmetry of the initial Hamiltonian 
with respect to px and py, see Eq. ([T]). In Figs. 2b, 
2c and 2d we show calculated contributions to the cur- 
rent for growing values of kox ■ For nonzero values of k^x , 
both jx and jy appear. It is seen that the frequency and 
the shape of ZB oscillations change with growing fcox val- 
ues. For growing kox different ?7^'^ become large and in 



t(fs) 



FIG. 2: Contribution of the Ki point of tlie Brillouin zone 
to the electric current in graphene at B=10T versus time, 
calculated for a Gaussian wave packet of the width dx = dy = 
L = 81.13A and various values of fco^;: a) kox = 0, b) kox = 
0.02A-\ c) kox = 0.035A-\ d) kox = 0.06k-\ Thick lines 
- jy{t), thin lines - jx{t). 



consequence different frequencies cj^ and oj^ dominate 
in sums pT]) and p^ . The striking feature seen in Figs. 
2a, 2b and 2c is, that the ZB is manifested by several 
frequencies simultaneously. This is a consequence of the 
fact that, as follows from Eq. in graphene the en- 
ergy distances between the Landau levels diminish with 
n, which results in different values of frequencies u>'^ and 
u!^ for different n. Thus, it is the presence of an ex- 
ternal quantizing magnetic field that introduces various 
frequencies into ZB. 

For sufficiently large values of kox only one frequency 
prevails, as shown in Fig. 2d. This is related to the fact 
that, as seen in the inset of Fig. 8 in Appendix B, the 
coefficients U^'^ in this regime have a pronounced max- 
imum around a specific value Umax- The dominant fre- 
quency is LUjnax = ^{.^Uraax + 1 - ^U^ax), which is sim- 
ply the cyclotron frequency for Umax, see Fig. 1. Thus, 
it might appear that the current shown in Fig. 2d corre- 
sponds to the simple classical cyclotron motion and the 
trembling motion is manifested only by the damping in 
time (see Appendix D). This is however not the case. 

In Fig. 3 we show the calculated current for kox = 
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FIG. 3: The same as Fig. 2d, but for larger time scale. FIG. 4: The same as in Fig. 2, calculated for fixed values 
Results for the first picosecond coincide with those shown in of kox ~ 0.035A~^ and widths dx = dy — Sl.lSA but at 
Fig. 2d. different magnetic fields. Results 4b are the same as those 

shown in Fig. 2d. 



0.06A-1 (the same as in Fig. 2d), but in a much larger 
time scale. It turns out that, after the ZB oscillations 
seemingly die out, they actually reappear. 

Thus, for all k^x values (including /cqk — the ZB 
oscillations have a permanent character, that is they do 
not disappear in time. This feature is due to the discrete 
character of the electron spectrum caused by a magnetic 
field. The above property is in sharp contrast to the 
no- field cases considered until present, in which the spec- 
trum is not quantized and the ZB of a wave packet has a 
transient character, see mathematical terms, 

due to the discrete character of the spectrum, averages 
of operator quantities taken over a wave packet are sums 
and not integrals, see Eqs. ^T7]i and (|18p. The sums 
do not obey the Riemann-Lebesgues theorem for inte- 
grals which guaranteed the damping of ZB in time for 
a continuous spectrum (see Ref. ^IS])- We consider the 
demonstration of a permanent character of ZB oscilla- 
tions for a discrete spectrum to be the main result of our 
present work. 

In Fig. 4 we show the ZB oscillations of the current 
for a constant wave vector /cqi: — 0.035A~^ at different 
magnetic fields. It can be seen that the intensity of a 
magnetic field has a dramatic effect on ZB: not only its 
frequency is changed but also the character of oscilla- 
tions. Lower magnetic fields are equivalent to higher kox 
values, since both lead to higher Landau levels involved 
(see Figs. 2d and 4a). Inversely, higher magnetic fields 
and lower kox values lead to lower Landau levels involved 
(see Figs. 2b and 4c). At very small magnetic fields there 
exist three regimes of ZB oscillations: the 'initial' oscil- 



lations dying out during a few femtoseconds (as shown in 
Fig. 2 of Ref. [l3|)j the second range of oscillations dy- 
ing out during several picoseconds, as shown in Fig. 2d, 
and the third range of permanent, somewhat irregular 
oscillations shown for t > Ips in Fig. 3. 

Finally, we calculate the displacements x{t) and y{t) 
of the wave packet. To this end we integrate expressions 
(fT7| and (|18p with respect to time using the initial con- 
ditions xq = x{0) and yo — y{0)- The results are plotted 
in Fig. 5 in the form of x — y trajectories for different 
initial wave vectors ko^. The direction of movement is 
clockwise. The trajectories span early times (Ips) after 
the creation of a wave packet. As mentioned above, the 
ZB oscillations do not die out in time which is reflected 
by infinite trajectories. In Fig. 5 the trajectories are 
shown around the point xq — yo ^ 0, whereas in reality 
the y component of the center is almost a linear function 
of kox : yo « koxL'^- 

All in all, the presence of a quantizing magnetic field 
has the following important effects on the trembling mo- 
tion. (1) For B ^ the ZB oscillation are permanent, 
while for B = they are transient. (2) For B many 
ZB frequencies appear, whereas for B = only one ZB 
frequency exists. (3) For _B ^ both interband and in- 
traband (cyclotron) frequencies appear in ZB; for B = 
there are no intraband frequencies. (4) Magnetic field 
intensity changes not only the ZB frequencies but the 
entire character of ZB spectrum. 

The results shown in Fig. 2 were obtained using the 
simplifying assumption about packet's width: dx = dy = 
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FIG. 5: Zitterbewegung trajectories of electron at the K\ 
point of the BriUouin zone during the first picosecond for var- 
ious values of fco^: a) k(,x = O.OlA"^, b) kox = 0.02A~^ c), 
= 0.035A-1 d), = 0.06A-\ 



FIG. 6: Intensity spectrum versus frequency during the first 
20 ps of motion of an electron described by a Gaussian wave 
packet having kox ~ O.OSSA"^ at B = 10 T. The frequencies 
and u'^ are defined in Eqs. (|17|l and (|18|l . 



L = [h/eBf/'^. This allowed us to use formula for 
the calculation of U^'^. However, in Fig. 4 we show 
the results obtained for constant dx — dy and variable 
B, for which we had to use general formula In 
all the calculations involving magnetic field, precise nu- 
merical values of the Hermite polynomials are required, 
see Eqs. and (P^ . For the results shown above we 
used the values of the first 400 Hermite polynomials, and 
checked their high precision using sum rule (jAl[) . We 
also considered the case of an initial electron momentum 
directed not along the x direction (as shown above), but 
also along the y direction. In this case the UZ'i coeffi- 
cients are imaginary, so that only the last two terms of 
Eqs. p?)) and (|18p come into play. The final results are 
similar but not identical to those quoted above. The rea- 
son is that, as already mentioned, the initial Hamiltonian 
(fT|) is not symmetric in and py momenta. When using 
Gaussian wave packet (I2ip . we assumed equal upper and 
lower components a„ = a; = 1/ \/2. This is in contrast 
to previous pa per s which usually took a„ = 1, a; = 
0) fill [T^ - [la . llSj l . This choice is somewhat arbitrary, 
it is determined by an experimental wave packet usually 
prepared by optical methods. One should keep in mind 
that the relative final amplitudes of Vx (t) and Vy (t) (and 
the resulting currents) depend on this choice via U^'jl co- 
efficients, see Eqs. (P5|) and If one chooses a„ = 1, 
a; = 0, the resulting motion coming from the Ki point 
of BZ is only along the y direction. 

It is of interest that phenomena analogous to those de- 
scribed above for electrons, occur also for photons. In 



particular, Hamiltonian ^ is similar to that describ- 
ing an interaction of atoms with electromagnetic radi- 
ation according to the so called Jaynes-Cummings model 
[stI . [38| . In particular, the collapse and revival of elec- 
tron's ZB oscillations, as illustrated in Figs. [2}d and[3l is 
predicted by the Jaynes-Cummings model for the num- 
ber of emitted photons ant is was observed in one-electron 
masers [391. 



VI. DIPOLE RADIATION DUE TO ZB 

Experimental possibilities of observing the trembling 
motion were considered previously 0, S [13 and we do 
not consider this problem again. The results shown in 
Figs. 2, 3 and 4 describe the electric current, which is an 
observable quantity. One could also try to see directly 
the displacement of charge, cf. [13, On the other 
hand, we argue below that the ZB should be accompa- 
nied by electromagnetic dipole radiation emitted by the 
trembling electrons. The oscillations r{t), as shown in 
Fig. 5, are related to the dipole moment — ef(t), which 
couples to the electroma gnet ic radiation. We shall treat 



the radiation classically i.e. we take the radiated 



transverse electric field to be [43 



er{t) sm{0) 
47reoc2 R ' 



(27) 



where eg is the vacuum permittivity, 9 is an angle between 
the direction of motion of a wave packet and the position 
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of the observer R. Integrating 5^ over the angle 9 one 
obtains the total radiated power given by the Larmor 
formula 



P = 



e2a2(t)2 



(28) 



where a is the acceleration averaged over the packet. To 
find P we calculate the acceleration from Eqs. ([T7| and 
by taking the time derivatives. The spectrum of the 
emitted radiation is obtained by the Fourier transform of 
the electric field. We have 



where 



1 X ^ ^ / imTt\ 
-ao+ 2^ a„cos I — ^ 1 

m=0 ^ ^ 



/ rmrt 



(29) 



= lim 



bm = lim 

T^oc 



cos 



sm 



imrt\ -> , 
h:^y^ir,t)dt, (30) 



and Oq = 0. For the numerical calculations we take a 
large period T = 20 ps. The intensity spectrum of oscil- 
lations is 



/(w,„) oc ^(C + ^) 



(31) 



0.4- 



0.2- 



0.0- 



B=10T 




0.02 



0.04 



0.06 (^-1) 0.08 



FIG. 7: Intensities of emission lines corresponding to uj^ 
(dashed line) , loq = ^ (dotted line) and to uj^ (solid line) ver- 
sus the wave vector kox of Gaussian wave packet at _B = 10 
T. 



The plot of I{o-!m) is given in Fig. [51 The strongest fre- 
quency peak corresponds to oscillations with the basic 
frequency lu = fl. The peaks on the high frequency side 
correspond to the interband excitations and are charac- 
teristic of ZB. The peaks on the lower frequency side 
correspond to the intraband (cyclotron) excitations. At 
higher magnetic fields we may expect smaller number 
of distinct frequencies, while for lower fields the classi- 
cal radiation will evolve toward a quasi-continuous spec- 
trum. In absence of Zitterbewegung the emission spec- 
trum would contain only the intraband (cyclotron) fre- 
quencies (see Appendix D). Thus the interband frequen- 
cies Lu^ shown in Fig. [5] are a direct signature of the 
trembling motion. It can be seen that the peaks are 
not drastically weaker than the central peak at lu ~ 
which means that there exists a reasonable chance to ob- 
serve them. 

In Fig. [7] we plot the dependence of the emitted power 
intensity on the initial wave vector kox for three lines (cal- 
culated for B = 10 T): the basic cyclotron line at w = 17, 
the intraband frequency ^3, and the interband frequency 
ujg . It is seen that the intensity of various emission lines 
depends differently on ko^. At small fco^^ values the basic 
line uj — Q dominates, but at fcox — 0.04A~^ the intensi- 
ties of various lines become comparable. The character- 
istic two maxima of ujq occur also for the other interband 
frequencies. We believe that the fcox dependence of the 
line intensities, as shown in Fig. [3 can serve as a signa- 
ture of ZB. 



The properties shown in Figs. [2] -[7] have been calcu- 
lated for the Ki point of the Brillouin zone in graphene. 
The main features applying to the K2 point of BZ are de- 
scribed in Appendix E. The above calculations are some- 
what idealized since they do not take into account the 
position of the Fermi energy in a given sample. Thus 
they correspond approximately to a situation with the 
Fermi energy relatively low in the valence band. Clearly, 
the frequencies corresponding to transitions with the fi- 
nal states below the Fermi energy are not possible. 

In general terms, the excitation of the system we pro- 
pose is due to the nonzero momentum hkox given to the 
electron. It can be provided by accelerating the electron 
in the band by or by exciting the electron with a nonzero 
momentum by light from the valence band to the con- 
duction band. The electron can emit light because the 
Gaussian wave packet is not an eigenstate of the system 
described by Hamiltonian ((!]). The energy of the emit- 
ted light is provided by the initial kinetic energy related 
to the momentum hkox- Once this energy is completely 
used the emission will cease. If the electron is described 
by a non-gaussian wave packet, all our results are quanti- 
tatively valid, only the intensity spectrum will differ from 
that shown in Fig. [51 We emphasize that the sustained 
character of ZB oscillations caused by a discrete energy 
spectrum makes graphene in a magnetic field probably 
the most favorable system for an experimental observa- 
tion of the trembling motion considered until present. 
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VII. SUMMARY 

We described the Zitterbewegung of electrons in solids 
in the presence of a magnetic field assuming that the 
electrons are represented by Gaussian wave packets. The 
system under consideration is monolayer graphene. ft is 
shown that the presence of a quantizing magnetic field 
has a profound influence on the Zitterbewegung. In par- 
ticular, the discrete energy spectrum in a magnetic field 
causes the Zitterbewegung to be sustained in time while 
for i? = the ZB has a transient character. In addi- 
tion, at i? ^ many ZB frequencies appear whereas at 
-6 = one deals with only one ZB frequency. For a given 
value of initial electron momentum, the magnetic field 
intensity affects not only ZB frequencies but the entire 
shape of the ZB spectrum. We consider and describe 
an electromagnetic radiation emitted by the trembling 
electrons. Graphene in a magnetic field seems to be a 
very favorable system for an experimental observation of 
Zitterbewegung . 

APPENDIX A 

We consider a sum rule for the coefficients U^'^ of Eq. 
([20]) . Let us calculate 1 = (/|/), 

n 71. s 

= 2j2{u::: + uk'^j. (At) 

n=0 

Factor 2 appears due to the summation over s. It is to 
be reminded that a'^ + af = f . The above sum rule can 
be used for a verification of the numerical values of J7"'". 



APPENDIX B 




FIG. 8: Contribution of the A'l point to the electron velocity 
Vy{t) versus time, calculated for the indicated parameters: (a) 
using full formula (|17|l (see Fig. 2d), (b) using the first inte- 
gral of the Poisson formula (|Bf |l . The two curves practically 
coincide. Inset shows coefficients C/n,n±i ~ Un±i,n ~ U{x) 
for kox = 0.12A~^. The frequency of oscillations corresponds 
to uJmax for Xmax = 46, See text. 

calculations of Eq. (H?]) for fcox = 0.f2A"^ The two 
curves practically coincide, apart from the small contri- 
butions of higher ZB frequencies present in exact formula 
()17p . The effective frequency of the motion is given by 
^max — ^(V Xmax ~\~ f \/ Xjnax)i where Xmax — 46 Cor- 
responds to the maximum of U{x), see inset. For larger 
times the second term in Eq. (jBf p is not negligible and 
full formula (|Bip is equivalent to Eq. P7)l . 



We calculate Vy{t) from Eq. (fT7| for a situation when 
U^'^ have a maximum for a large value of 71. In our 

model there is C/"_'^ = aaapUn,n- We use the Poisson 
summation formula for an estimation of the velocity av- 
erage disregarding and ap coefficients. Upon replacing 
Um,n by a continuous variable U{x) and approximating 
Un,n±i = Un±i.n ~ U {x) , the term with V~ in Eq. (fT7|) 
vanishes. Then 

/•oo 

Jo 

GO ^ 

-|-4u^ / sm{ujxt)U{x) cos{2TTxl)dx, (Bf) 
1=1 



where ujx = ^{s/x -f I — \/x). For sufhciently small times 
we may disregard the second term and Vy{t) is given 
by the first integral in Eq. (jBf p . In Fig. [5] we show 
the results of the integration compared with the exact 



APPENDIX C 

Here we consider briefly the gauge aspects. According 
to the general theory [44] , if one introduces a new gauge 
by means of an arbitrary function A(r), the new vector 
potential \s A' ^ A + VA, and the new scalar potential 
is A'q ~ Aq + dA/dt. Then the wave function changes as 
_ g(ie/?i)A^^ ^YiQ gauge invariance for an operator 
O — 0{A, Aq) means 



{^\OiA,Ao)\^) 
This leads to 



{^'\d{A',A'„m 



diA,Ao) 



,d{A,Ao) 



(CI) 



(C2) 



If, instead of the gauge A = {—By, 0,0), we take A' = 
(0,Bx,0), which gives the same magnetic field, we have 
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A = xy/L"^, so that *'(a;,y) = e'^2'/^''^'(x, y). Using 
prescription (|C2p and calculating 



e"^/^ ,Py 



(C3) 



one shows that relation (|C1|) is satisfied also for Hamil- 
tonian ([I]) in the new gauge. 



Thus an electron represented by a Gaussian wave packet 
having the initial momentum hko^ moves on a circular 
orbit with the cyclotron frequency lUc without attenua- 
tion. A similar result is known for a one-dimensional 
wave packet moving in a parabolic potential. On the 
other hand, the motion illustrated in Fig. 2d is damped 
during the first picosecond which is an another manifes- 
tation of Zitterbewegung. 

APPENDIX E 



APPENDIX D 

We consider here the motion of a wave packet in the 
presence of a magnetic field according to the Schrodinger 
equation. For 2D Hamiltonian there is |n) = \nk.j;) = 
e''=-"H„(0e-«'/V(V2^C„), and £;„ = hujdn + 1/2), 
where ujc = eB/m. The velocity average Vx{t) is 

Ut) = - V (/|n'fc^>(nfc,|/)e'-=("'-")*(n'fc;|7r,|nfc,), 

(Dl) 



m 

n,n'k^ ,k' 



and similarly for Vy(t). Since 



{nk'Jn^lnkx) 
we have 



V2L 



{\^Sn',n-l + VrH-iSn' ,n+l) , 

(D2) 



Here we consider contributions related to ZB of elec- 
trons at the inequivalent point K2 of the Brillouin zone. 
The form of Hamiltonian at the K2 point is somewhat 
controversial, various authors give different expressions. 
According to Refs. [4^ and [4y] the Hamiltonian H' is 



H' 





-TTrr + ITT, 



y 







(El) 



i.e. H' = —H^. The eigenvectors of H' are 



\nkxs)' = 



74^ V 1) 



(E2) 



v,{t) = --=— / dkxY.^f\n-l){n\f)^e-'^-^' + 
V2Lm J 

h I" °° 
- / dfc.V(/|n + l)(n|/)^Ame'-=^ (D3) 



i.e. they differ from those given by Eq. ([5]). The quantum 
velocity dH' /dpx — —uux and dH' /dpy ~ +uay. Thus, 
the X component of the velocity changes sign, while the 
y component remains unchanged. Repeating the calcu- 
lations we obtain again Eqs. ([TT]) and (HH]), in which the 
coefficients are 



There is ^/n\n — 1) = d\n) and \/n + l\n + 1) = a+jn), 
and we calculate 



V2Lr 



dkx{{fa\f)e-'-^' + {fa+\f)e'-^'), 

(D4) 

and similarly for Vy (t) . For a one-component wave packet 
of Eq. pT|) the integrals indicated in Eq. (jD4| can be 
done analytically. We finally obtain 



= 

^ n 


n.n- 


- f/"'" m jj'-''- 


-1 + '^n-l.n' 


f± ^ 




hi '^n+l,n "T '^n,n 


_|_ ttI,1 
-1 '^n-l.n' 


4± - 


^ n.n 


1 TtI,u _|_ ttU,1 

1 '^n.n ^ "-^n+l.n-l 


+ "-^n-l^n-hl' 


Bt = 


_JJU,1 

7i,n 


'^n,n "T '-'ri+l,n-l 


+ '^n-l^ri+l- 



Vx {t) 

Vy{t) 



2m 
.fiko. 



hki 



Ox 



COs{uJct), 



2 m 



hki 



Ox 



(E3) 



If we are interested in the electric current, we should 
add the velocities of the two inequivalent points of the 
BZ. As a consequence, the x component of the velocity 
vanishes while the y component nearly doubles. 
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